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We review how a high-statistics observation of the neutrino signal from a future galactic core- 
collapse supernova (SN) may be used to discriminate between different neutrino mixing scenarios. 
Most SN neutrinos are emitted in the accretion and cooling phase, during which the flavor-dependent 
differences of the emitted neutrino spectra are small and rather uncertain. Therefore the discrimina- 
tion between neutrino mixing scenarios using these neutrinos should rely on observables independent 
of the SN neutrino spectra. We discuss two complementary methods that allow for the positive iden- 
tification of the mass hierarchy without knowledge of the emitted neutrino fluxes, provided that the 
f3-mixing angle is large, sin 2 813 3> 1CP 5 . These two approaches are the observation of modulations 
in the neutrino spectra by Earth matter effects or by the passage of shock waves through the SN 
envelope. If the value of the 13-mixing angle is unknown, using additionally the information encoded 
in the prompt neutronization v e burst — a robust feature found in all modern SN simulations — can 
be sufficient to fix both the neutrino hierarchy and to decide whether #13 is "small" or "large." 
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I. INTRODUCTION 



Despite the enormous progress of neutrino physics 
in the last decade, many open questions remain to be 
solved. Among them are two, the mass hierarchy — 
normal versus inverted mass spectrum — and the value 
of the 13-mixing angle #13, where the observation of 
neutrinos from a core-collapse supernova (SN) could 
provide important clues [H, 0, @] ■ Schematically, the 
neutrino emission by a SN can be divided into four 
stages: Infall phase, neutronization burst, accretion, 
and Kelvin-Hclmholtz cooling phase. During the infall 
phase and neutronization burst only f e 's are emitted, 
while the bulk of neutrino emission is released in all 
flavors in the two latest phases. Whereas the neu- 
trino emission characteristics of the two initial stages 
are basically independent of the features of the pro- 
genitor, such as its mass or equation of state (EoS) of 
the core, the details of the neutrino fluxes during the 
accretion and cooling phases may significantly change 
for different SN models. Therefore, a straightforward 
extraction of oscillation parameters from the bulk of 
the SN neutrino signal seems hopeless. Only features 
in the detected neutrino spectra that are independent 
of unknown SN parameters should be used in such an 
analysis. 

In this talk we discuss the potential of the two most 
promising sources for such features: the modulations 
in the P e spectra caused by the Earth matter 0, d, II] 
or by the passage of shock waves through the SN enve- 
lope [EUTil- Moreover, we show how the detection of 
the neutronization v e burst could help to break possi- 
ble degeneracies in the case that #13 is still unknown. 
In Tab. U it is shown the dependence of these three 
observables on three different neutrino mass schemes. 



TABLE I: The presence of Earth-matter and shock wave 
effects in the v e spectra and of the v e burst for different 
neutrino mixing scenarios. 



Case Hierarchy sin 2 #13 Earth Shock v e burst 



A 


Normal 


> 


10" 3 


Yes 


No 


No 


B 


Inverted 


> 


1(T 3 


No 


Yes 


Yes 


C 


Any 


< 


10- 5 


Yes 


No 


Yes 



II. IDENTIFYING SIGNATURES OF THE SN 
SHOCK WAVE PROPAGATION 

The neutrino spectra F Vi arriving at the Earth are 
determined by the primary neutrino spectra F%. as 
well as the neutrino mixing scenario, 

F 1/i (E,t) = YlPd E ,t) F 2 j ( E ,t), (!) 
j 

where pji is the conversion probability of a i/j into Vi 
after propagation through the SN mantle. The prob- 
abilities pji are basically determined by the number 
of resonances that the neutrinos traverse and their 
adiabaticity [2l|. Both are directly connected to the 
neutrino mixing scheme. In contrast to the solar case, 
SN neutrinos must pass through two resonance lay- 
ers: the H-resonance layer at pn ~ 10 3 g/cm 3 cor- 
responding to Am5; tm , and the L-resonance layer at 
Pl ~ 10 g/cm 3 corresponding to Am|. Whereas the 
L-resonance is always adiabatic and in the neutrino 
channel, the adiabaticity of the H-resonance depends 
on the value of #13, and the resonance shows up in 
the neutrino or antineutrino channel for a normal or 
inverted mass hierarchy respectively [l[ . 

During approximately the first two seconds after 
core bounce, the neutrino survival probabilities are 
constant in time and in energy for all three cases A, 
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FIG. 1: Density profile at several instances after core 
bounce. The resonance layers pu and pl are also shown §. 




FIG. 2: Antineutrino survival probability pv c v e (E,t) as 
function of energy at different times averaging in energies 
with the energy-resolution of Super-Kamiokande. A den- 
sity profile with a forward and reverse shock have been 
used, from Ref. 



B, and C. However, at t s» 2 s the H-resonance layer 
is reached by the outgoing shock wave, see Fig. [TJ 
The way the shock wave passage affects the neutrino 
propagation strongly depends on the neutrino mixing 
scenario: cases A and C will not show any evidence 
of shock wave propagation in the observed v e spec- 
trum, either because there is no resonance in the an- 
tineutrino channel as in scenario A, or because the 
resonance is always strongly non-adiabatic as in sce- 
nario C. However, in scenario B, the sudden change 
in density breaks the adiabaticity of the resonance, 
what leads to a time and energy dependence of the 
antineutrino survival probability pa s v c (E, t) , see FigJ3J 
The presence of the shocks results in the appearance 
of bumps, one in the case of a unique forward shock 
or two if an additional reverse shock is present, in 
Pv c v c (E, t) at those energies for which the resonance 
region is passed by the shock waves. All these struc- 



tures move in time towards higher energies, as the 
shock waves reach regions with lower density, leading 
to observable consequences in the v e spectrum. 

A useful observable to detect effects of the shock 
propagation is the average of the measured positron 
energies, (E e ), produced in inverse beta decays D e + 
p — > n + e + , the most important neutrino signal ex- 
pected in a SN. In Figs. [3] and HI we show (E e ) together 
with the one sigma errors expected for a megaton wa- 
ter Cherenkov detector and a SN in 10 kpc distance, 
with a time binning of 0.5 s. For the neutrino fluxes 
we use the parametrization suggested in Refs. [HJdll], 



F° Vj {E) 



$o(i/j) (a,- + 

{EoM) r( aj + i) 

( - E 



E 



exp 



(Eaivj)) 



(2) 



The values of the parameters used in the two models, 
Gl in Fig. and G2 in Fig. gj are given in Tab. HU 
The averaged energy is assumed to decrease linearly 
after 5 s. 



TABLE II: The parameters of the used primary neutrino 
spectra models motivated from the SN simulations of the 
Garching (G1,G2) and the Livermore (L) group. The aver- 
aged energy is given in MeV and the pinching parameters 
a Pe and a Sx take the values 4 and 3, respectively. 

Model (E {P e )} (E (u x )) go(gg)/go(i^) 



L 


15 


24 


1.6 


Gl 


15 


18 


0.8 


G2 


15 


15 


0.5 



The effects of the shock wave propagation are 
clearly visible, and are independent of the assump- 
tions about the initial neutrino spectra @. Moreover, 
it is not only possible to detect the shock wave prop- 
agation in general, but also to identify the specific 
imprints of the forward and reverse shock versus the 
forward shock only case. The signature of the reverse 
shock is its double-dip structure compared to the one- 
dip of a forward shock only. The observation of the 
details of the shock wave structure will be, though, 
tightly related to the absence of strong density fluctu- 
ations after the shock front na, UM . 



III. EARTH-MATTER EFFECTS 

As it has been previously mentioned, before the 
shock wave reaches the H-resonance layer the depen- 
dence of the conversion probability in the cases A, 
B and C on the neutrino energy E is very weak. 
However, if neutrinos cross the Earth before reaching 
the detector, the conversion probabilities may become 
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FIG. 3: The average energy of vp — > ne + events binned 
in time for a static density profile, a profile with only a 
forward shock, and with forward and reverse shock. Gl 
model was assumed for the neutrino fluxes. The error bars 
represent 1 a errors in any bin, from Ref. 0. 
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FIG. 4: Same as Fig. |3] for the G2 model. From Ref. 



energy-dependent and induce modulations in the neu- 
trino energy spectrum. These modulations may be 
observed in the form of local peaks and valleys in the 
spectrum of the event rate <jFP plotted as a function 
of 1/E. These modulations arise in the antineutrino 
channel only in cases A and C. Therefore its obser- 
vation would exclude case B. This distortion in the 
spectra could be measured by comparing the neutrino 
signal at two or more different detectors such that the 
neutrinos travel different distances through the Earth 
before reaching them 0, [TtJ • However these Earth 
matter effects can be also identified in a single detec- 
tor 0,i]. 

The net v e flux at the detector may be written in 
the form 

F? = sin 2 9 12 Fl + cos 2 12 Fl 

7 

+ AF°J2 A i^ 2 (kiV/2), (3) 

i=l 



where y is the "inverse energy" parameter y = 
12.5 MeV/E, AF° = (F^ - F D o), and A l depend 
only on the mixing parameters. 




FIG. 5: Averaged power spectra in the case of a large scin- 
tillator detector for different SN models, Gl (top panels) 
and L (bottom panels), and distances travelled through 
the Earth, from Ref. [g]. 

The last term in Eq. Q is the Earth oscillation term 
that contains up to seven analytically known frequen- 
cies ki in y, the coefficients AF°Ai being relatively 
slowly varying functions of y. The first two terms in 
Eq. ^ are also slowly varying functions of y, and 
hence contain frequencies in y that are much smaller 
than the fcj. The frequencies ki are completely inde- 
pendent of the primary neutrino spectra, and can be 
determined to a good accuracy from the knowledge 
of the solar oscillation parameters, the Earth matter 
density, and the position of the SN in the sky Q . The 
latter can be determined with sufficient precision even 
if the SN is optically obscured using the pointin g ca - 
pability of water Cherenkov neutrino detectors [la ]. 
The power spectrum of N detected neutrino events is 



G(k) 



1 

N 



N 



e iky ' 



(4) 



In the absence of Earth effect modulations, G(k) has 
an average value of one for k > 40. The region k < 40 
is dominated by the "0-peak", which is a manifesta- 
tion of the low frequency terms in Eq. ([3]). Identi- 
fying Earth effects is equivalent to observing excess 
power in G(k) around the known frequencies fcj. In 
Fig. [5]we show the averaged power spectra in the case 
of a 32 kton scintillator detector for neutrinos travel- 
ling only through the Earth mantle (left panels) and 
traversing both mantle and core (right panels). It is 
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possible to observe how in the former case only one 
peak is present, whereas in the case that neutrinos go 
through the core more frequencies arise. The model 
independence of the peak positions may be confirmed 
by comparing the top and bottom panels. 

Since in the real world the presence of fluctuations 
in the signal, see Fig. [SJ will spoil any naive theoreti- 
cal peak, we need to introduce a prescription to carry 
out the analysis. One possibility is to consider the 
area around the expected position of the peak. The 
area under the power spectrum between two fixed fre- 
quencies A; m in and fc max is on an average (fc m ax — fcmin)- 
In the absence of Earth effects, this area will have a 
distribution centered around this mean. The Earth ef- 
fect peaks tend to increase this area. The confidence 
level of peak identification, p a , may then be defined as 
the fraction of the area of the background distribution 
that is less than the actual area measured. 




V 



FIG. 7: Comparison of pgs as a function of nadir an- 
gle n for a 32 kton scintillator (SC) and a megaton water 
Cherenkov (HK) detector 0. 
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FIG. 6: Realistic power spectrum from a single simula- 
tion, from Ref. [(| . 



In Fig. [3 we compare P95 obtained with a 32 kton 
scintillator detector and a megaton water Cherenkov 
detector assuming a SN at 10 kpc and the Gl model 
for neutrino fluxes. In the latter case, as neutrinos 
travel more and more distance in the mantle the peak 
moves to higher k values, and due to the high k 
suppression, the efficiency of peak identification de- 
creases. When the neutrinos start traversing the core, 
additional low k peaks are generated and the efficiency 
increases again. The identification of Earth matter ef- 
fects excludes case B, and is thus complementary to 
the observation of shock wave effects. 



IV. NEUTRONIZATION v e BURST 

If the value of #13 is unknown, a degeneracy ex- 
ists between case A and C. Both scenarios predict the 
same v e signature in a water Cherenkov detector, and 
therefore the previous two observables are not useful 
to disentangle them. In this case, the additional in- 



formation encoded in the v e neutrinos emitted during 
the neutronization burst can fix the range of On as 
well as the neutrino mass hierarchy. 

The prompt neutronization burst takes place during 
the first ~ 25 ms after the core bounce with a typical 
full width half maximum of 5-7 ms and a peak lumi- 
nosity of 3.3-3.5 xlO 53 erg s _1 . The striking similar- 
ity of the neutrino emission characteristics despite the 
variability in the properties of the pre-collapse cores 
is caused by a regulation mechanism between electron 
number fraction and target abundances for electron 
capture which establishes similar electron fractions in 
the inner core during collapse. This leads to a con- 
vergence of the structure of the central part of the 
collapsing cores and only small differences in the evo- 
lution of different progenitors until shock breakout. 
The small dependence of the neutronization burst on, 
e.g., the p roge nitor mass can be verified in Fig. [8] (cf. 
also Refs. pj|20|). 

Taking into account that the SN will be likely ob- 
scured by dust and a good estimation of the distance 
will not be possible, the time structure of the detected 
neutrino signal should be used as signature for the 
neutronization burst (l9| . Since the event number in 
current and proposed charged-current v e detectors is 
not high enough to allow for a detailed time analysis, 
we discuss only the case of a megaton water Cherenkov 
detector. Here one has to consider the v e elastic scat- 
tering on electrons, which is affected by several back- 
grounds like inverse beta decay or reactions on oxygen. 
This background can be substantially reduced by us- 
ing angular and energy cuts, as well as Gadolinium to 
tag neutrons from inverse beta decays. The sample of 
elastic scattering events still contains the irreducible 
background of scattering on electrons of other neu- 
trino flavors than v e , but this contamination does not 
affect the possibility to disentangle the different neu- 
trino scenarios (l9j . 
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The time evolution of the signal depends strongly 
on the neutrino mixing scheme. In case A, the v e 
survival probability is close to zero, and therefore the 
peak structure observed in the initial v e luminosity 
is absent. On the contrary, in case C, 30% of the 
original v e remain as v e whereas 70% are converted 
into v x . Since the cross section of v e on electrons is 
much larger than that of v x , the signal is dominated 
by the contribution of v e . These i/ e 's follow the time 
evolution of L„ c , and thus lead to a clear peak in the 
signal. 



In Fig. [9] we show the expected neutrino signal from 
t = —5 to 18 ms for different progenitor masses, and 
for the mixing scenarios A and C. The peak structure 
can be clearly seen in case C, but not in case A fl9| . 
Including recent improvements of the electron capture 
rates or uncertainties in the nuclear equation of state 
has only little effect on the neutronization peak com- 
pared to the size of the statistical fluctuations. There- 
fore the observation of a peak in the first milliseconds 
of the neutrino signal would rule out case A. 




tpb (ms) 

FIG. 8: Neutrino luminosities as function of time for dif- 
ferent progenitor masses, from Ref. 



V. SUMMARY 

A reliable determination of neutrino parameters us- 
ing SN neutrinos should be independent from the pri- 
mary neutrino fluxes produced during the accretion 
and cooling phase of the SN. Earth-matter effects and 
the passage of SN shocks through the H-resonance 
both introduce unique modulations in the neutrino 
energy spectrum that allow one their identification 
without knowledge of the primary neutrino spectra. 
While the observation of Earth-matter effects in the 
v e energy spectrum rules out case B, modulations in 
the v e time spectrum identify case B. If the value of #13 
would be known to be large, then the neutrino mass 
hierarchy would be identified. Otherwise, the detec- 
tion of the neutronization v e peak — a robust feature of 
all modern SN simulations — can break the remaining 
degeneracy between A and C. 
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FIG. 9: Number of ve — > ve events per time bin in a 
megaton water Cherenkov detector for a SN at 10 kpc for 
cases A (dashed lines) and C (solid lines) and for different 
progenitor masses. Statistical errors are also shown for the 
15 M case [ijj. 
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